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ABSTRACT 



Aims. We attempt to determine the optical depth of spiral galaxy disks by a statistical study of new Tully-Fisher data from the ongoing 
KLUN+ survey, and to clarify the difference between the true and apparent behavior of optical depth. 

Methods. By utilizing so-called normalized distances, a subsample of the data is identified to be as free from selection effects as 
possible. For these galaxies, a set of apparent quantities are calculated for face-on positions using the Tully-Fisher diameter and 
magnitude relations. These values are compared with direct observations to determine the mean value of the parameter C describing 
the optical depth. 

Results. The present study suggests that spiral galaxy disks are relatively optically thin tb ~ 0.1, at least in the outermost regions, 
while they appear in general to be optically thick tb > 1 when the apparent magnitude and average surface brightness are studied 
statistically. 

Key words, galaxies: spiral - galaxies: photometry - galaxies: ISM - galaxies: fundamental parameters 
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1. Introduction 

The optical depth of spiral galaxy disks has been debated exten- 
sively with widely diverse results. The optical depth is an im- 
portant parameter to determine because it affects how galaxies 
appear at different viewing angles. By understanding the incli- 
nation dependence of total magnitude, surface brightness, and 
diameter, inclination corrections for these observed quantities 
can be determined, and spiral galaxies inclined at different an- 
gles and located at different distances can be related. The opti- 
cal depth of spiral galaxies is at particular interest in supernova 
surveys because of t he ne ed to correct supernova magnitudes 
dHatano et all 1 19971 Il998l) and supernova rate determinations 
kiello & Patatt l2005t iBotticella et all 120081) for the effects of 
extinction. Internal extinction is another important parameter re- 
lated to selection bias in extragalact ic Cepheid distance determi- 
nations (Paturel & Teerikorpii 120051) . 

The classical approach is to derive overall disk extinction 
from the effects of disk inclination on observables such as mag- 
nitude, diameter, and average surface brightness. If there is no 
dust in the spiral galaxy disks or the dust extinction is negligibly 
small, the apparent magnitude of the disk appears to be constant, 
regardless of viewing angle because no material absorbs the ra- 
diation and the galactic environment is so diffuse that individual 
star images, such as point-like sources, do not overlap each other. 
However, the apparent surface brightness increases as the incli- 
nation increases because the area of the ellipse projected by the 
galaxy image decreases. For the optically thick case, the bright- 
ness quantities vary in an opposite way. If there is a lot of dust 
in the spiral galaxy only the outer layer of the stars can be seen 
or at least it dominates the brightness observations. The depth of 



this star layer remains constant inspite of the inclination angle 
and therefore so does the surface brightness. However, as the in- 
clination increases and the apparent surface area decreases, the 
apparent total brightness will decrease also. 

Several statistical studies applying the classical or simi- 
lar methods have claimed that spiral galaxy disks are opti- 
cally t h ick at least in their in ner parts. These include Valentiin 
(1990), Burst ein et alj d 199 lh including distance information, 
Giovanelli et al. (1994) studying surface bri ghtness shifts, 
Peleti er et al. (1994) using multiple wavebands, iBottinelli et all 
(119951) utilizing Tully- Fisher relation (here after TF relation) with 
normalized distances, iTullv et ail d!998l) w ith multiple wave- 
bands and volume limited galaxy g roups, Shao et all (T2007) 
studying the luminosity function, and lDriver et alj ( 20071) using 
a inclination relation model for the galaxy components. 

Although the studies cited above and related papers im- 
ply that spiral galaxy disks should be considered to be opti- 
cally thick, conflicting results have been obtained from sta tisti- 
cal studies. From diameters, Huizinga & van Albadi (119921) in- 
ferred th at galaxy disks are optically neither completely thick 
nor thin. iDavies et alj ([1993) concluded that it is impossible to 
determine opti cal depth w i th stat istical studies because of se- 
lection effects. I Jones et ail d 19961) criticized the way in which 
galaxy samples had been used to obtain a volume-limited statis- 
tical sample. 

Statistical studies of the optical depth of spiral galaxies with 
samples of several hundreds or thousands of galaxies involve 
multiple selection effects, therefore inspecting only a few galax- 
ies closely, represents a completely different way to approach 
the problem. Investigations limited to edge-on galaxies usually 
led to the conclusion that the galaxies are optically thin. These 
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included the studies of Bosn ia et al, 



d 19921) andlBvunl 



1991. 



which compared rotation curves, and lXilouris etail d 1 997l fl998l 
1999), which applied various exponential fits to observations of 
stars and dust. 

For face-on spirals, a more direct method has been used: one 
observes background objects through the target galaxy. These 
studies suggest that the galaxies are optically rather thin, with 
optically thicker spiral arms and thinner interarms. A simpler 
fo rm of this method is to observe eithe r gal axy pairs, as achieve d 
bv lAndredakis & van der Kruitl (1 19921) andlBerlind et al.l(ll997l) . 
or multiple pai r s, as observed by IWhite et alj ( 20001) and 
iDomingue et al.l d2000t) . Another method of observing back- 
ground objects through spiral ga laxies is the so-called syn- 
thetic field method introduced by Gonzal ez et al.l (JT998), who 
studied the amount of distant galaxies observable through the 
parts of spiral galaxy disks compared with reference fields. 
iHolwerda et al.l d200 5a) improved the method w ith automation, 
and a st udy of several galaxies was presented bv lHolwerda et al] 
(I2005bl) . 

In contrast, the Milky Way disk is consid ered to be optically 
thin. A t the Galactocentric solar distance, IKoppen & Vergelvl 
(1998) inferred a value Ay = 0.1 for the face-on ab- 
sorption. However, even at high latitudes there appears to 
be directions in which there is higher local extinction e.g. 
iBerdvugin & Teerikorpil(l2002l) . 

In this study, we aim to determine the optical depth of the 
spiral galaxy disks and inclination corrections for diameter and 
magnitude in the new TF data sample obtained from the still on- 
going KLUN+ (Kinematics of the Local Universal survey. The 
structure of the paper is the follo wing. In Sect. 2, the method 
used for the old KLUN sample by Bottinelli litail d 19951) is pre- 
sented. In Sect. 3, the method is applied to a new sample and the 
results are given. In Sect. 4, an additional test of the method is 
presented using a large galaxy sample obtained from the extra- 
galactic database HyperLeda^. In Sect. 5, the difference between 
true and apparent optical depth is discussed and in Sect. 6, a 
summary is given. 



2. TF relation in the classical approach 

2.1. Influence of selection effects 

Almost every observed galaxy sample is limited by either di- 
ameter or magnitude or both, due to the limitations of the ob- 
servational program and the telescope adopted. Therefore, it is 
useful to summarize briefly how the classical approach, utilizing 
diameters and magnitudes versus inclination, is affected by the 
inevitable limits in angular size and/or apparent magnitudes. It 
is helpful to consider first the idealized case of a galaxy cluster 
consisting of a single type of galaxy, with a small scatter in linear 
diameter (or absolute magnitude). This can be easily generalized 
to field galaxies, which are localized into different bins of radial 
velocity, and therefore different "clusters", at different distances. 
If the sample is diameter limited and galaxies are opaque, the 
isophotal diameter does not depend on inclination: the sample is 
therefore representative of all inclinations and the diameter limit 
does not affect the (zero) slope. Classes of smaller galaxies just 
remain unobserved. 

In non-opaque case, the diameter increases with inclination; 
the low inclination part of the sample is then less complete, 
which produces too shallow a slope in the logD versus B/A 




B/A 

Fig. 1. Samples of spiral galaxies from the HyperLeda database 
separated into two groups with apparent magnitude 13.0 < m < 
13.2 above and 16.5 < m < 16.7 below. If limited by angular 
diameter, i.e. line logics = 0.8, the sample would lack face-on 
(axis ratio B/A — > 1) galaxies in the non-opaque case, and the 
slope would be artificially too shallow. 




1 http://klun.obs-nancay.fr/KLUN+/pagel .html 

2 http://leda.univ-lyonl.fr/ 



B/A 

Fig. 2. Samples of spiral galaxies from the HyperLeda database 
separated into two groups with apparent diameter log D25 > 0.90 
and logics < 0.70. If limited with diameter, i.e. line logZ^s = 
0.8, the sample would be complete in a non-opaque case. 



diagram, and therefore too high an opacity. The axis relation 
R = A/B is used if logarithmic scaling is needed where B is the 
minor axis and A the major axis of the galaxy. The sample com- 
pleteness example constructed from the HyperLeda database is 
shown in Fig.Q] 

In a magnitude-limited sample, the selection effect is oppo- 
site. In the opaque case, the magnitude depends on inclination 
(edge-on galaxies are fainter), so the high inclination part of 
the sample is less complete, which again produces too shallow 
slope in the m versus B/A diagram, and in this case artificially 
low opaqueness. In fully non-opaque cases, all inclinations are 
equally completely represented (see for comparison Fig. [2] for 
the trivial case). 

To avoid these effects, one may use the TF relation for mag- 
nitude or diameter. One constructs the Hubble parameter H ver- 
sus normalized kinematic distance diagram to find the unbi- 
ased plateau, where the sample is unaffected by the diameter (or 
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magnitude) limi t and s hould be representative of all inclinations 
dBottinelliet all 119951) . 



2.2. The relevant quantities 

For every spiral galaxy the apparent Hubble constant is calcu- 
lated by an infall model around the Virgo cluster, for which 



Ho 



(VoVir + Vo)4ta 



(1) 



where the distance r originates in M — m - 25 - 5 log r, in which 
the absolute magnitude M is evaluated by the TF relation 



M = fl(log V m - 2.2) + b , 



(2) 



where a and b are the TF parameters. Alternatively, one may also 
use the diameter TF relation. In what follows, the parameters of 
the diameter relation are denoted by a' and b' . 

A short review of the equation s derived and the gala xy model 
applied to the spiral galaxies by Bo ttinelli et al.l d 1995b is given. 
The inclination dependence of apparent magnitude given in the 
form m - mo + A(R) is 



= m Q - 2.5 log()fc + (1 - k)R 2c(1+0 



2/K)-l 



)■ 



(3) 



We ad opted the constants K — 0.1 from Fouqu e~& Paturell 
( 1985 b, when the morphologica l type w as T > 2, and k — 0.16 
from Sim ien & de Vaucouleursl (|l986) as a mean value when 
2 < T < 7. When the apparent surface brightness is given by 



/z = m + 5 log L>25 , 



(4) 



the inclination dependence of apparent surface brightness is eas- 
ily derived in the form [i-[io + B{R) as 



(X=(i - 2.5 log(/UT 2c + (1 - k)R' 



? (0AC/K)-1 



(5) 



2.3. Results bv lBottinelli et al.\ \1995) 

The initial sample of Bottinel fi et al.l d 1995b included 4535 spiral 
galaxies. Using the diameter TF relation, an unbiased subsample 
was identified. This subsample consisted of 332 galaxies. The 
value obtained for the subsample C = 0.07 + 0.05 corresponds 
to t S3 0.9 (see Fig. |3]for the dependence of r on C). However, 
it should be noted that the optical depth of the disk is expected 
to decrease as the radius increases. Therefore, parts of the inner 
disk may have r > 1 . The limitation of the plateau subsample did 
not change the result significantly. For the unlimited "biased" 
sample, the result C = 0.05 + 0.02 was obtained. 

3. New sample 

3.1. Basic sample 

The HI data used in this study were acquired as part of the KLUN 
project, based on observations with the Nangay radio telescope 
during the years 1990-2006, for the purpose of distance and pe- 
culiar velocity measurements. Approximately 5000 HI profile^] 
have been gathered at Nancay and corresponding data published 



Data tables and Hl-profiles and corresponding com- 
ments are available in electronic form at the CDS via 
anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/Abstract.html They are also provided 
through Simple Spectra Access and partly available at NED 
(http://nedwww.ipac.caltech.edu/!. 




Fig. 3. Dependence of optical depth t on the parameter C used 
in our method. 



in the KLUN series paper s dTheureau et all Il998l I2005L l2007i 
llMIPatarel et alll2003bb . 

All KLUN HI spectra acquired by the Nangay radio tele- 
scope antenna were reviewed and assigned a quality code ac- 
cording to the shape of their 21 -cm line profile. This study 
dGuillardl 12004 and references above) has enabled us to iden- 
tify efficiently a large number of TF outliers due to morpholog- 
ical type mismatch or Hi-confusion in the elongated beam of 
Nangay. 

The photometric input catalogue was compiled from the 
Hyperleda extragalactic database: 

- Accurate coordinates (of typical accuracy lower in value 
than 2 arcsec), diameter, axis ratio, and position angle were 
extracted from the analysis of the Digitalized Sky Survey 
(DSS); Diameters and axis ratios were corrected to the RC2 
system at the limiting surface brightness of 25 mag arcsec -2 
dPaturel et all l2003al) . 

- B-band magnitudes were extracted both from a compila- 
tion of Z?-band photometry available in the litera ture and, in 
a hom ogeneous way, directly from the DSS dPaturel et all 
2000). Each magnitude is the result of a weighted mean 
over all references, after correction to the RC3 system (see 
iPaturel et all 1 19971) . 

- To complete the Gal actic extinction corre ction, a g , we used 
values adopted from Schlege l et al.l d 1998b . 

The morphological classification in HyperLeda is based 
on visual inspection. As for the other parameters, each 
type is obtained as a weighted mean of the T parameter 
dVaucouleurs et all 1199 ll) over existing references in the litera- 
ture. This number T is between -5 and -2 for ellipticals, between 
-2 and 1 forlenticulars, 1 to 8 for spiral galaxies, and 9 and 10 for 
irregular. The DSS image of each KLUN target was inspected by 
eye before HI observation, to confirm its spiral galaxy nature and 
determine whether a bar was present. The Sa type spiral galax- 
ies were excluded from the sample because of the possibility that 
the bulge of the galaxy may dominate the light detected during 
observations. Galaxies close to the Galactic plane \b\ < 15 ° with 
large galactic extinction were excluded. 

The adopted method of normalized distances utilizes the 
Virgo cluster infall model, with parameters RAy u = 12.5 h, 
#vir = 12°, and Vovh = 980km/s. For the infall speed of 
the Local group, a v alue of vo = 150km/s is used as in 
Botti nelli et al.l d 19951) . Galaxies at small angular distance from 
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Fig. 4. Inclination dependence of magnitude for the initial sam- 
ple of 772 galaxies with the mean trend plotted. 
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Fig. 5. Inclination dependence of diameter for the initial sample 
of 772 galaxies with the mean trend plotted. 

the Virgo cluster 9 < 30 ° and Virgo members are excluded so 
that unambiguous values of d^ n = r /ryj r can be determine d from 
the kinematic velocity field model (Bo ttinelli et al.L[T995f) . Here 
r is the distance of the sample galaxy and ry\ r the distance of the 
Virgo cluster so c4in and later on d n are dimensionless quantities. 
All of these restrictions further limit the size of the sample to 772 
individual galaxies. From here on, this is referred to as the initial 
sample unless otherwise specified. The inclination dependence 
of magnitude and diameter for the initial sample is presented in 
Fig.[4]and[5] respectively. 

In Fig. |H it can be seen that the magnitude distribution of 
the brightest and closest galaxies is rather smooth. Regardless of 
how the distant galaxies are selected, the distribution of galax- 
ies with the lowest m v alues in the sample differs s ignificantly 
e.g., from the ESO-LV dLauberts & Valentiinl fl989) sample in 
which the brightest edge-on galaxies are on average 1 — 1.5 mag 
fainter than the brightest face-o n galaxies. This can be s een in 
sev eral opacity studies s uch as [Huizinga & van Albadal (l992) 
and lBurstein et alj dl99ll) . 

The initial sample can be taken to be complete for a diameter 
limit of log D25 = 1.1 and magnitude limit of m-\ « 15. The sam- 
ple therefore contains more distant galaxies than the data used by 



Botti nelli et al.l ( 1 1995b . The sample is not limited by magnitude, 
although, galaxies with logics < 1.1 are restricted from the 
sample. The diameter limit is also used in Eq. ©. From the di- 
ameter distribution in Fig. [31 it can be seen that t he sample lacks 
faint e dge-on galaxies. This was also noted by Bott inelli et al.l 
(119951) . who attributed the lack to a systematic error in which 
the inclinations of small galaxies were measured to be too smal l. 
The same problem was discussed by iGiovanelli et al.l ( 1 19941) . 
who claimed that it was a problem caused by the resolution 
limit of observations. Wh atever the source of the problem is, we 
use the same limitations as lBottinelli et alJ (TT995). who excluded 
galaxies with log(7?) > 0.8 {B/A < 0.16). Face-on galaxies with 
\og(R) < 0.07 (B/A > 0.85) are excluded because they lack ac- 
curately measured rotation speed and are of little use to the TF 
relation. 

This limits the sample to 495 galaxies. This is referred to as 
the basic sample unless otherwise specified. For the basic sam- 
ple, initial TF relation parameters were determined by fitting 
slopes to a relation between linear diameter and the maximum 
rotational velocity with a value of a' = 1.1, and to a relation 
between absolute magnitude and rotational velocity with values 
a = -6.3 and b = -19.52. A zero point V for the diameter is 
insignificant in the adopted method because the face-on diame- 
ters determined by the TF relation are scaled to the observations 
with a fitting constant as seen in Eq. (0. 

3.2. The diameter limited subsample 

The basis of the normalized distance method is to extract at each 
distance, a subsample free of selection effects, i.e. a subsam- 
ple for which the average TF distance is its true mean distance. 
Normalized distances for the basic sam ple are derived from th e 
kinematic distances by an equation from lBottinelli et all (1 19951) 

d B = d km lO a ' {2J - losV ™ } • 10- clo ° (A/B) • 10 0094fl -« • 10 logD '- L3 , (6) 

where logDi is the diameter limit for the sample determined 
above and a g is the galactic extinction. The optical depth pa- 
rameter has an initial value of C = 0. The average value of 
the Hubble constant is calculated for galaxies at the same nor- 
malized distance d n . As expected from the selection effect that 
removes absolutely small galaxies, the value of Ho is seen to in- 
crease at large normalized distances. This effect is seen in Fig. [6] 
at d n > 6. With smaller values, the Hubble constant has an 
average value of Ho » 60kms~ 1 Mpc~ 1 forming an unbiased 
sample consisting of 122 spiral galaxies. This is referred to as 
the diameter-limited (DL) subsample unless otherwise specified. 
The diameter TF parameter with a value a' — 1 .2 was determined 
for the DL subsample. 

For the DL subsample, the apparent reference diameters that 
the galaxies should have if seen face-on are determined with the 
TF relation 

log D = - log d kin + a' log Vm + est . (7) 

The parameter C that provides the average optical depth of the 
DL subsample is determined from the equation 

logD 25 =logD + ClogOR). (8) 

The inclination dependence of log(D25/£>o) is presented in 
Fig. [7] The line fit to the data is completed with a least squares 
method and has a slope value of C = 0.24 + 0.05. It should be 
mentioned that for the method used here, of all observed values 
the accuracy of the apparent diameter has the most significant 
effect on the result. 



E. Kankare et al.: Optical depth of spiral galaxies 



x°8 




Fig. 6. Hubble constant versus the normalized distance for the 
basic sample. The curve represents the constant value of the 
Hubble parameter for the subsample when d n < 6 and the in- 
crease of the parameter with a linear fit when d n > 6. 
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Fig. 8. Inclination dependence of magnitude difference m — mo. 
The curves from top to bottom repre sent the function A(R) with 
C = 0.01 for these data, C = 0.04 by lBottinelli et al.l(ll995l) . and 
C = 0.24 for the DL subsample, respectively. 
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Fig. 7. log(D25/Do) versus \og(R) for the DL subsample. The 
upper line is a fit for the data with a slope C = 0.2 4. The lower 
line represents the result of iBottinelli et alJ (1 19951) with a slope 
C = 0.07. 



The value obtained for parameter C corresponds in terms of 
optical depth to r ~ 0. 1 which would imply that spiral galaxies 
are extremely optically thin. However, as expressed in Sect. 2.2, 
it is expected that the optical depth decreases as the radius in- 
creases. The result C = 0.24 for the DL subsample is true only 
for both the diameter and the outer rim of the disk, and does not 
provide any information about the optical depth in the inner parts 
of the disk or even about the average optical depth for the entire 
disk. This is because the method relies on diameter observations 
that are affected mostly by the optical depth in the outer parts of 
the galaxy disk. 

3.3. The magnitude limited subsample 

The behavior of the brightness quantities and derived parameter 
value C can be tested with the functions presented in Eqs. (f3]l and 
(0. Since the equations adopt brightness observations for the 
entire galaxy disk projected on the sky, there is a need for another 



subsample. The initial sample of 772 galaxies is now limited as 
before but by the magnitude completeness limit instead of the 
diameter limit providing another basic sample of 574 galaxies 
with TF parameter values of a = -5.7 and b = -19.67. The 
normalized distance is now determined for magnitude purposes 
by equation 



dn = 4inl0 a2a(1 °8 V »- Z7) • 10' 



0.2A(«) 



10 0.2a g . 10 



-0.2(m T -13) 



(9) 



where inclination correction has an initial value of A(R) = 
2.5 log/?, which corresponds to C = and an opaque disk. The 
behavior of the Hubble constant versus the normalized distance 
is similar regardless of the way in which the distance is de- 
termined from the kinematic distance, and again galaxies with 
the normalized distance d n > 6 are excluded providing another 
subsample of 124 galaxies with essential parameter values of 
a = -6.5 and b = -19.19. For each galaxy, the apparent refer- 
ence magnitude if seen face-on is calculated using the TF rela- 
tion 



mo = 5 log (iki n + a log V m + est . 



(10) 



The function A(R) and the ML subsample are presented in Fig. [8] 
The best-fit function to the data provides as low a value as C = 
0.01 + 0.05, corresponding to t > 1, which agrees f airly well 
with the result of C = 0.04 from lBottinelli et al.l dl995l) . 

Similarly the function B(R) is presented in Fig. [9] The ap- 
parent face-on surface brightness fio that the galaxy should have 
is derived from Eqs. (|4]i and ( fTQb . The best-fit description of the 
data is C = 0.02 + 0.07 (corresponding again to t > 1). Since the 
total brightness difference between face-on and edge-on galaxies 
is large, it dominates the surface brightness difference regard- 
less of the increase in apparent diameter due to the optically thin 
outer galaxy disk. 

3.4. Other subsamples 

Additional information can be obtained by restricting the sam- 
ple in different ways. The optical depth of the basic sample can 
be determined by applying limits to the diameter as achieved 
for the DL subsample above. The method pr ovides a value 
C = 0.07 + 0.03, which agrees with the results of Bott inelli et al.l 
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Fig. 9. Inclination dependence of surface brightness difference 
fi - fiQ. The curves represent the function B(R). The upper curve 
is a fit for the data with a value C = 0.02. The lower curve 
represents the result from the DL subsample with a value C = 
0.24. 
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Fig. 10. log(D25/Z)o) versus log(7?) for the basic sample. The 
upper line represents the result from the DL subsample with a 
slope C = 0.24. The lower line is a fit for the data with a slope 
C = 0.07. 



d 1995b . However, the difference between the DL subsample and 
the basic sample for the new data is not unexpected. The result 
derived for the DL subsample implied that spiral galaxy disks 
are optically thin at least in terms of the inclination dependence 
of the diameter; this would imply that the basic sample limited 
to certain values of diameter is incomplete and lacks faint, low- 
inclination-angle, galaxies, as stated in Sect. 2.1. A graphical 

presentation is s een in Fig. [TUl 

We note that Botti nelli et alj (1 19951) considered galaxies with 
HI profile observations that had only one peak or a low S/N ra- 
tio. To avoid a possible source of error only galaxies with high 
quality HI profile data were included in the samples used above 
to determine the optical depth. For comparison the data samples 
could be expanded to include also low-quality observations. This 
increases the diameter-limited basic sample to 835 galaxies with 
TF parameters of a' = 0.9, a = -4.7, and b - -19.67. When 
limited by the normalized distance, a subsample of 332 galaxies 
is obtained with a TF parameter of a' = 1.1. The method then 
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Fig. 11. log(D25/A>) versus \og{R) for the DL subsample in- 
cluding low-quality data. The lines from top to bottom represent 
C = 0.24 for th e basic DL subsample, C = 0.18 for these data, 
and C = 0.07 by lBottinelli et alJ (Il995l) . respectively. 
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Fig. 12. log(D25/A)) versus \og(R) for the DL subsample with 
different morphological types separated as Sab-Sb (o) galaxies, 
Sbc-Sc (•) galaxies, and Scd-Sd (A) galaxies. See text for details. 



provides a slope value C = 0.18 + 0.03 for the opacity parameter 
with the fit presented in Fig. [TT] The result corresponding to a 
value t w 0.2 is relatively close to the value obtained for the DL 
subsample; it therefore appears unable to explain the di fferences 
between our results and those of Botti nelli et alj d 19951) . 

Finally, the DL subsample, as can be seen in Fig. Q~2l is di- 
vided into closer morphological types. The dispersions and re- 
sults for the categories are: 46 type 2 < T < 3.5 Sab-Sb galaxies 
with a value C = 0.24 + 0.10, 51 type 3.5 < T < 5.5 Sbc-Sc 
galaxies with a value C = 0.23 ± 0.08, and 25 type 5.5 < T < 7 
with a value C = 0.25 + 0.07. There appears to be no signifi- 
cant differences between morphological classes, all types being 
optically very thin at least in their outer parts. 



4. HyperLeda sample 

Because the data sample used in Sect. 3 was small, a more sta- 
tistical study was completed for a far larger galaxy sample that 
included HI observations obtained from the HyperLeda. No re- 
strictions were placed on the quality of the radio observations. 
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0.2 



0.4 



0.6 



0.8 



log(R) 



Fig. 13. log(D25/Do) vs. \og(R) for the HyperLeda sample. The 
upper line is a fit for the dat a with a slope C = 0.1 8. The lower 
line represents the result of Botti nelli et alj (I 1995b with a slope 
C = 0.07. 



The sample includes 9580 spiral galaxies from which galaxies 
with a small angular distance to the galactic plane and the Virgo 
cluster were excluded, such that the sample was similar to the 
initial sample in Sect. 3.1. Using similar limitations, a selection- 
effect-free subsample of 989 galaxies with d n < 5 was obtained 
with the TF relation. With the method of normalized distances 
andEqs. 0}-© and ©-I©, a slope value of C = 0.18+0.03 was 
obtained, which was equivalent to r » 0.2 for the log(D2s/A)) 
versus log/? relation (see Fig. [TSl . This result correlates well 
with the results for the subsample in Sect. 3.2. 

For a magnitude-limited subsample, a sample of 345 spiral 
galaxies was obtained from the HyperLeda data with the normal- 
ized distance d a < 1.5. The method represented by Eqs. ©-(O 
and d9T>-(TT0T> corresponds to a parameter value C = 0.02 + 0.05 
from the total magnitude relation and C = 0.00 + 0.06 from the 
surface brightness relation with both corresponding to r > 1. 
This again agrees well with the results in Sect. 3.2. All results 
obtained from different samples for the optical depth parame- 
ter C with error limits, used TF parameters, normalized distance 
limits, and sample sizes are presented in Table Q] 

5. Discussion 

5.1. Effects that increase apparent optical depth 

It appears essential to divide the true optical depths of spiral 
galaxy disks from those values inferred observationally from sta- 
tistical arguments. As individual low-inclination-angle galaxies 
have been studied extensively utilizing different techniques, such 
as observing galaxy pairs or using the synthetic field method, 
it appears that face-on galaxy disks are quite optically thin. 
However, statistical results infer the presence of thicker disks. 
There are at least three possible reasons for this disagreement be- 
cause the statistical studies concentrate usually on the inclination 
dependence of apparent brightness quantities and the brightness 
difference between face-on and edge-on galaxies in a classical 
way: light scattered by dust, the geometric distribution of the 
absorbing dust, and the uncertainty in diameter. 

Light scatteri n g can affect the apparent optical depth. 
iBvun et all d 1994b. IPi Bartolomeo et all (Il995l) . IXilouris et all 
1997b . and iBaes & Deionghd (12001 ) studied several models of 



light scattering in galaxy disks. Their models indicated that even 



if spiral galaxies are relatively optically thin, part of the light 
emitted by the galaxy itself scatters from the thin dust and will 
increase the apparent luminosity of the low inclination galax- 
ies. When individual galaxies are studied, this has no effect, for 
example, because background objects are observed through the 
galaxy under study. The scattering from the light emitted from 
the background obj ects can also be ne glec ted as demonstrate d 
by studies such as lBerlind et al.l(ll997l) and lWhite et alj (|2000). 
However, if optical depth in galaxies is studied statistically, in- 
ternal scattering will cause apparent brightness to have an in- 
clination dependence, which may be interpreted classically as 
evidence for optically thick galaxy disks. 

The second effect is the geometric distribution of dust. Close 
studies (see Sect. 1) concluded that spiral galaxies are thin 
not only over the entire disk, but are divided into optically 
thicker spiral arms and thinner interarms. This type of absorb- 
ing medium structure implies that a classical sandwich model is 
too simplistic for describing spiral galaxies. This is, because the 
ratio between the projected areas of the thicker spiral arms and 
the galaxy disk increases as the inclination angle increases. This 
will again increase the apparent brightness of face-on galaxies 
relatively and cause disks to behave as if they were optically 
thicker if studied statistically. However, it should be noted that 
this effect is minimized by the fact that optically thicker spiral 
arms are also the brightest parts of the disks and therefore the 
dominant components in any surface or total brightness studies. 

The third possible effect for poor quality data is the error 
in measured diameters. Giovanell i et al.l (1 19941) compared scale 
lengths and surface brightness diameters of observed galaxies 
with standard table values concluding a decrease in the observ- 
able surface brightness levels. Therefore, visual determination of 
the apparent diameter without using the apparent surface bright- 
ness levels would underestimate the size (and also therefore the 
apparent total brightness) of the edge-on galaxies. 

These effects increase the apparent total brightness differ- 
ence between face-on and edge-on galaxies and therefore have 
an impact on the apparent optical depth when inclination de- 
pendence of total brightness is studied statistically. The impor- 
tance of this division into real and apparent optical depth comes 
from the need for accurate inclination corrections in observed 
galaxies to determine the values that these would have if seen 
face-on. This is because close optical depth studies of individual 
galaxies do not necessarily provide accurate inclination correc- 
tion factors. Interestingly, the conclusion that galaxies are really 
optically thin but behave as if they were thicker c ontradicts the 
conclu sions for the sandwich model presented by iDisnev et al.l 
(1989), which claimed that spiral galaxies could be optically 
thick but behave as if they were thin because of the dominance 
of the outer crust of the stars. This underlines the fact that the 
statistical approach depends on the model required to interpret 
the distributions of and correlations between observable quan- 
tities. Furthermore, it can be unable to detect any spherically 
symmetric dust distribution. For example, th ere is evidence fo r 
reddening in the halo components of galaxies dTeerikorpii 12002 ) . 

5.2. The value of C in statistical studies 

The results of our statistical study are in agreement, indirectly, 
with the above effects, although they cannot explaine entirely 
the high values of optical depth. The results also agree with the 
widely accepted idea that the outer regions of the disks are opti- 
cally thinner than the inner parts. Thirdly, the inclination correc- 
tions for apparent magnitude and appare nt surface brightness d o 
not differ significantly from the results of Bottine lTi et al.l (1 1995b . 
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Table 1. Limiting normalized distances, number of galaxies, 
used TF parameters, and results for different samples. 



N Name < rf n Size a' a b C C ar r 



1 


DL subsample 


6.0 


122 


1.2 


-7.5 


-18.92 


0.24 


0.05 


0.1 


2 


ML subsample (m) 


6.0 


124 




-6.5 


-19.19 


0.01 


0.05 


> 1 


3 


ML subsample (ji) 


6.0 


124 




-6.5 


-19.19 


0.02 


0.07 


> 1 


4 


DL basic sample 




495 


1.1 


-6.3 


-19.52 


0.07 


0.03 


0.9 


5 


DL subsample" 


6.0 


332 


1.1 


-6.5 


-19.12 


0.18 


0.03 


0.2 


6 


DL subsample b 


5.0 


989 


1.0 


-5.7 


-19.53 


0.18 


0.03 


0.2 


7 


ML subsample b (m) 


1.5 


345 




-5.1 


-19.13 


0.02 


0.05 


> 1 


8 


ML subsample b (p) 


1.5 


345 




-5.1 


-19.13 


0.00 


0.06 


> 1 



a Same as in N = 1, now including low-quality data. 
b HI data from the HyperLeda. 



A difference with earlier work is appearent in inclination cor- 
rections for the isophotal diameter, which could be resolved by 
more accurate observations and more pre cise TF relation param- 
eters. Using the same parameter values as lBottineili et all ( 1 19951) . 
the same method (see Sect. 2 & 3) yields a subsample consist- 
ing almost entirely of the galaxies in the DL subsample in Sect. 
3.2 and corresponding to a slope value of C ~ 0.20. This has no 
significant effect on the value of t but decreases the difference 
between the slope values. 

It should also be noted that the value derived from the 
diameter correctio n agrees well with the v alue C = 0.235 
used i n the RC2 ( Vaucouleurs et all Il976l) . whereas the re- 
sult of IBottinelli et all (1 19951) agrees fairly well with the RC3 
dVaucouleurs et all 1 199 ll) value of C — 0.0 for spiral galaxies. 
In view of the present results, we recommend that one checks 
whether the choice between these values of C influences the fi- 
nal results. 



5.3. Implications on high-z studies 

Local TF studies are being conducted increasingly more us- 
ing infrared data. The progre ss of more local TF research has 
been towards infrared bands. | Masters et ail d2006h studied the 
SFI++ sample for / and lMasters et al.l d2003l 120081) the 2MASS 
Extended Source Catalog for J, H, and K bands. Ther e have 
also been TF surveys at higher reds hifts dWeiner et all [2006; 
IChiu et all l2007t iKassin et all 120071) . These studies aim to de- 
termine the evolution in the luminosities and kinematics of spiral 
galaxies. Our results provide a starting point for applying ex- 
tinction corrections to these data. However, the characteristics 
of these samples differ considerably from ours in terms of wave- 
length and selection criteria. Furthermore, the isophotal radii of 
the more distant galaxies only reach their inner regions, because 
of the dimming of the surface brightness caused by the redshift. 
The derived inclination angle is also expected to be less accurate, 
which does not affect the corrections of quantities with slight 
inclination dependence (diameter in our study), but does cause 
extra scatter in quantities with stronger inclination dependence 
(e.g. magnitude in our study). Obviously, a statistical study of a 
more distant sample would be a more challenging project, which 
may provide results that different from ours. 



6. Conclusions 

In a statistical study of opacity, we have analyzed data from the 
ongoing KLUN+ survey with a method that includes TF and 
kinematical distance information. 



1 . Our analysis has found that spiral galaxy disks are quite opti- 
cally thin, at least in their outer regions, with an optical depth 
value tb *0.1 obtained using t he diameter TF relation . This 
differs from the conclusion of IBottinelli et all (1 19951) . who 
inferred that the outer regions were optically thick. 

2. However, using apparent magnitude and apparent surface 
brightness as test quantities, it appears that the entire disk 
behaves statistically as a thick component with an apparent 
optical de pth value Tb > 1. Thi s is in approximate agree- 
ment with IBottinelli et all d 1995b and with the magnitude- 
inclination correction adopted by HyperLeda. 

3. On the other hand, various studies of face-on galaxy disks 
have found optically thin disks. This apparent contradiction 
with statistical studies based on inclination angle may be at 
least partly due to light scattering, geometric distribution of 
the dust, and incorrect measurements of the apparent diam- 
eter. These factors are expected to have an effect on the in- 
clination correction values, which can be resolved accurately 
only by larger, statistical methods. 
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